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- - -. ABSTRACT 
O 

o 

CN ' We report on the abundance of interstellar neutral nitrogen (NI) for 30 sight- 

lines, using data from the Far Ultraviolet Spectroscopic Explorer (FUSE) and the 
Hubble Space Telescope (HST). NI column densities are derived by measuring the 
■ equivalent widths of several ultraviolet absorption lines and subsequently fitting 

those to a curve of growth. We find a mean interstellar N/H of 51 ±4 ppm. This 
is below the mean found by Meyer et al. of 62+g ppm (adjusted for a difference 



Or 



in /-values). Our mean N/H is similar, however, to the (/-value adjusted) mean 



> 

of 51 ±3 ppm found by Knauth et al. for a larger sample of sightlines with larger 
hydrogen column densities comparable to those in this study. We discuss the 



question of whether or not nitrogen shows increased gas-phase depletion in lines 
of sight with column densities logA^(H to t) > 21, as claimed by Knauth et al.. 
The nitrogen abundance in the line of sight toward HD 152236 is particularly 
interesting. We derive very small N/H and N/O ratios for this line of sight that 



2 ■ may support a previous suggestion that members of the Sco OBI association 

^ ■ formed from an N-deficient region. 
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1. INTRODUCTION AND BACKGROUND 

Many atomic species in the interstellar medium (ISM) show enhanced depletions relative 
to hydrogen in dense, highly reddened environments. For the most part, however, neutral 
nitrogen (NI) has appeared to be an exception. Because nitrogen is not expected to be a 
significant source of material for dust grains (Snow & Witt 1996), and the abundances of 
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molecules such as N 2 also appear to be small in the ISM (Knauth et al. 2003), it is not 
surprising that nitrogen depletions may be small. Therefore, if enhanced nitrogen depletion 
in a given line of sight is detected, it would challenge current models of interstellar dust 
composition and/or nitrogen chemistry. 

In diffuse interstellar clouds, which are permeated by a radiation field whose high-energy 
cut-off is 13.60 eV, nitrogen is primarily in the form of NI because its ionization potential 
of 14.53 eV is slightly higher than HI. Hence, the analysis of nitrogen depletions depends 
on this assumption that nitrogen in HI regions remains largely unaffected by ionization and 
therefore the total nitrogen abundance relative to hydrogen, N/H, is accurately approximated 
by JV(NI) /[iV(HI) + 27V(H 2 )]. 

One of the earliest studies of interstellar NI, by Lugger et al. (1978), attempted to 
empirically determine the /-values of many NI absorption lines in the Copernicus wavelength 
range. The unknown /-values were determined from one line of sight where several lines 
with known oscillator strengths were on the linear portion of the curve of growth. This new 
information was used to determine column densities for eight other sightlines. Their study 
found varying depletions of nitrogen, but with large errors. The determined abundances 
were all within a factor of four of the best solar value at that time (Withbroe 1971). 

Ferlet (1981) observed many additional lines of sight with Copernicus. Using the /- 
values determined by Lugger et al., Ferlet found the nitrogen abundance to be independent 
of reddening (i.e. total hydrogen content) in the line of sight. Ferlet also found the average 
nitrogen abundance to be consistent with a solar abundance within the large uncertain- 
ties (the solar abundance has been subsequently revised; see discussion below). Also using 
Copernicus observations, Hibbert, Dufton, & Keenan (1985) found N/H to be mostly in- 
dependent of reddening. A small exception was that less reddened lines of sight appeared 
to have slightly more nitrogen depletion than reddened lines of sight. Hibbert et al. were 
unable to rule out the possibility that hydrogen column densities were overestimated in the 
unreddened cases. 

In further Copernicus observations, York et al. (1983) determined the average value 
of N/H in the ISM to be 50% of the solar value, and also noted it was independent of 
reddening. While the York et al. value for N/H in the ISM and the solar value of N/H could 
be reconciled somewhat within the uncertainties, their results were the first to suggest that 
nitrogen in the ISM must be somewhat depleted if the solar abundances of various elements 
are appropriate standards to apply to the ISM. It is worth noting that the lines of sight 
considered by York et al. had a visual extinction of 1 magnitude at the most. Additionally, 
York et al. used the /-values presented by Cowan et al. (1982), which are much larger than 
the /-values determined by Lugger et al. (1978). Therefore, the nitrogen column densities of 
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York et al. were systematically smaller than the nitrogen column densities of prior studies. 

Meyer, Cardelli, & Sofia (1997) used the Goddard High Resolution Spectrometer (GHRS) 
onboard the Hubble Space Telescope to study lines of sight with low visual extinction (A v < 
1 mag). Their study found the value of N/H in the ISM to be 80% of the Grevesse & Noels 
(1993) solar value (75 ±4 ppm for the ISM compared to 93 ± 16 ppm). However, within the 
uncertainties, the results of Meyer et al. could not rule out the possibility of either a solar 
N/H abundance or an abundance of N/H that is 2/3 solar. The case of a 2/3 solar N/H 
abundance would have been consistent with depletions observed at the time for O/H and 
Kr/H. However, the idea that many ISM abundances are 2/3 of the solar abundance is no 
longer accepted (Sofia & Meyer 2001), based largely on revisions to the solar abundances 
that have brought ISM abundances such as O/H into closer but not necessarily complete 
agreement with the solar values (see also the discussion and references in Jensen et al. 2005). 
NI column densities determined by Meyer et al. were based exclusively on measurements of 
the weak 1160 A doublet. Therefore, it is important to note that the uncertainties of the 
Meyer et al. study included these /-values of the weak doublet, determined by Hibbert et al. 
(1985). Adjusted for the more recent /-values calculated by Tachiev & Froese Fischer (2002) 
and used in this paper, the Meyer et al. mean is 62^ ppm (see §3.3 for more discussion of 
these /-values). For the sake of direct comparison, throughout the rest of this paper we 
quote this adjusted mean. 

Moos et al. (2002) used the Far Ultraviolet Spectroscopic Explorer (FUSE) to study the 
H, D, N, and O content of several nearby sightlines and discovered variations in the N/H 
and N/O abundances. Moos et al. concluded that this difference is due to local ionization. 

Recently, Knauth et al. (2003) measured NI column densities in lines of sight with 
Ay > 1 mag using data from FUSE and the Space Telescope Imaging Spectrometer (STIS) 
onboard HST. They found the first significant evidence of enhanced nitrogen depletion as 
total hydrogen column density increases. Specifically, Knauth et al. found nitrogen abun- 
dances less than the Meyer et al. mean with respect to both hydrogen and oxygen for lines 
of sight with logiV(H tot ) > 21. Knauth et al. also placed limits on the N 2 abundance that 
ruled it out as a significant sink of NI in the lines of sight with enhanced depletions. The 
weighted average N/H abundance Knauth et al. found, including lines of sight from Meyer et 
al. (with some final N/H ratios slightly adjusted), Jenkins et al. (1999), and Sonneborn et al. 
(2000) — the latter two sources using IMAPS and different absorption lines — is 59 ± 2 ppm; 
adjusting the Knauth et al. and Meyer et al. results for the Tachiev & Froese Fischer (2002) 
/-values, this mean abundance is 54 ± 2 ppm. The /-value adjusted weighted average of the 
17 lines of sight analyzed by Knauth et al. is 51 ± 3. As with the Meyer et al. mean, we 
quote this adjusted mean throughout the rest of this paper. 
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The purpose of this paper is to examine the NI content of many additional moderately 
reddened lines of sight, some of which have more H 2 and a higher A v than typical lines of 
sight in previous studies. 

2. OBSERVATIONS AND DATA REDUCTION 

Data were taken from the FUSE and HST archives. The sightlines were chosen from a 
sample observed by FUSE as part of a molecular hydrogen survey conducted by Rachford et 
al. (2002) and continued in B. L. Rachford et al. (in preparation). The original FUSE sample 
included 37 sightlines, and relevant HST data were available for eight of these sightlines. For 
three sightlines (HD 166734, HD 281159, and NGC 2264 67, i.e. Walker 67) we could not 
reliably measure any NI absorption lines. Additionally, we detected the 1134 A triplet of 
NI in the spectrum for the lines of sight toward HD 43384, HD 62542, HD 164740, and HD 
210121. However, the column densities that might be derived from profile fitting or a curve- 
of-growth method based on this moderately-strong triplet alone would be very uncertain 
because the profile is relatively insensitive to column density and 5-value. Systematic error 
might also be a concern because a variety of multiple-component velocity dispersions could 
conceivably replicate the profiles, but with a very different total column density than derived 
from a single-component profile fit or curve of growth. 

Thus, our sample is reduced to 30 lines of sight where we can calculate a column density 
for NI, including 26 lines of sight for which we can reliably calculate N/H based on direct 
measurements of both HI and H 2 . In the remaining four cases, direct measurements of H 2 
are available, but HI must be inferred. In three of these remaining four cases we use the 
relationship presented by Bohlin et al. (1978) to infer the total hydrogen column density 
from Eb~v'i m the fourth case we quote a literature value (Hanson et al. 1992) where V(HI) 
is inferred from E B -v (see §4.2). Basic stellar data for the stars of these sightlines are given 
in Table 1. Values for hydrogen column densities and reddening parameters along the lines 
of sight are given in Tables 2 and 3, respectively. 

2.1. FUSE Data 

Data are collected with the FUSE satellite by four channels, two with a lithium fluoride 
detector (LiFl and LiF2) and two with a silicon carbide detector (SiCl and SiC2). Each 
channel contains two data segments, designated as A and B, covering adjacent wavelength 
regions. The lithium fluoride channels cover the wavelength region from 989-1188 A, while 
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the silicon carbide channels cover the region from 905-1104 A. 

The data were reduced with the CALFUSE pipeline, version 2.4.0. Data from different 
data segments are not coadded; rather, the data segment with the best signal-to-noise (S/N) 
is selected, with other data segments providing a check for consistency. 

The resolution element of FUSEis ~ 15 — 20 km s _1 . In fitting the non-Gaussian profiles 
exhibited by many of the lines in this study, we used a Gaussian point-spread function 
(PSF) with full-width at half-maximum (FWHM) of 15kms _1 . Our primary objective is 
to determine the equivalent widths of the various absorption lines (as opposed to deriving 
column density or b- value from individual lines); therefore, the results of our fits do not 
depend strongly on the exact PSF that is used. 

2.2. HST Data 

Archived HST data taken by STIS were available for ten sightlines: HD 24534, HD 
27778, HD 37903, HD 147888, HD 185418, HD 192639, HD 206267, HD 207198, HD 207538, 
and HD 210839. All of these data sets provide information on the NI triplet at 1200 A except 
the data for HD 24534. 

We used on-the-fly calibrated STIS data (Micol et al. 1999). In many cases we had 
multiple observations. We coadded observations of similar echelle order, adding errors in 
quadrature. We used empirically measured STIS PSF's (S. V. Penton, private communi- 
cation) appropriate to the grating and aperture of the observations in the fitting of the 
non-Gaussian profiles that are always exhibited by the 1200 A triplet. The PSF's we use are 
definitively non-Gaussian and in other contexts might be worth more discussion; however, 
they are narrow enough compared to the profiles that the fits and resulting equivalent widths 
for these very strong lines do not depend strongly on the precise nature of the PSF that is 
used. 

For many of these sightlines, IUE data are available, covering a wavelength region similar 
to that of HST. However, the few lines that are in the IUE wavelength region are blueward 
of Lyman-a and subject to considerable uncertainty that we deem unacceptable for inclusion 
in this study. 
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3. OBSERVED NI ABSORPTION LINES 

We follow the same procedure outlined in our previous work on 01 (Jensen et al. 2005). 
Specifically, we measure the equivalent widths of as many neutral nitrogen absorption lines as 
possible, then fit the measured equivalent widths to a curve of growth represented by a single 
velocity dispersion (i.e. b- value). The equivalent widths of at least three absorption lines 
are required to produce a meaningful solution for the two free parameters of column density 
iV and 6-value (measuring two absorption lines produces a unique but poorly-constrained 
solution; see Jensen et al.). In particular, to better constrain the solution and check for 
possible systematic errors due to the single 6-value assumption, it is necessary to measure 
absorption lines with a range of /-values, at least some of which are on the linear portion of 
the curve of growth. 

The following is a brief description of the NI absorption lines that were used in our 
analysis. Wavelengths, /-values, damping constants, and the instrument that was used for 
the observations of each of these lines are summarized in Table 4. All atomic data in this 
section are taken from the compliation of Morton (2003); the /-values relevant to this study 
are originally found in Tachiev & Froese Fischer (2002). In one line of sight (HD 210839) 
we detected all of the absorption lines used in this study; as an example, these absorption 
lines are shown in Figure 1. 

3.1. Lines Between 950 A and 960 A 

There are three multiplets between 950 A and 960 A that are observable. One is a 
doublet at 951 A. The stronger component of this doublet was detected and measured in 
nine lines of sight, and the weaker component was found in five of those nine lines of sight. 
Another doublet is found at 959 A. The stronger component of this doublet is found in 13 
lines of sight, while the weaker component is not detected in any lines of sight. These three 
lines are all relatively weak and resolved velocity structure is not apparent; therefore, the 
broad instrumental profile of FUSE dominates the observed profile, and Gaussian fits were 
appropriate in all cases (even if the intrinsic profile is slightly saturated). 

In between these two doublets is a moderately strong triplet found at 953 A. In our 
lines of sight, all of which possess large amounts of H2, the reddest component of the triplet 
is masked by H2 absorption. However, the other two components are observed. Either 
Gaussian or Voigt profiles were used to fit the lines, as appropriate to the observed profiles 
of the lines. 
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3.2. 1134 A Triplet 

The triplet found at 1134 A is the most commonly observed set of absorption lines in 
this study, as the lines are the strongest NI transitions found in FUSE spectra and they also 
lie in a spectral region with relatively high S/N. These lines were detected in every sightline 
for which we were able to measure any NI absorption lines. There are several complications 
in analyzing these lines. First, in some spectra there are apparent stellar features that must 
be fit to get an accurate representation of the continuum. Secondly, the separation of the two 
weaker (and bluer) components of the triplet is so small that the lines blend in many cases. 
This problem is enhanced in cases of high b- value and/or column density. Thirdly, there 
are cases where the lines are clearly saturated but the cores do not reach zero. This implies 
the presence of scattered light, contamination from telluric emission lines, or inaccuracies 
in the CALFUSE pipeline reduction. There are also known instrumental features (caused 
by exposure to telluric emission) in this location of the LiF2A spectra that have increased 
over time; thus we rely on LiFIB data, except in the earliest spectra where both LiFIB and 
LiF2A provide good data. 

We account for the non-zero cores by fitting a quadratic through the cores of the three 
components of the triplet, and subtracting this from the spectrum. This is a valid method 
if the non-zero cores are caused by a smooth, continuous source of contamination — which 
may or may not be true of scattered light or a problem in the reduction pipeline. Telluric 
emission would complicate matters, because while it is fairly broad (~ 0.3 A for each line), 
it would not be as uniform as our correction assumes. However, we are able to constrain the 
expected level of telluric NI emission by observing other strong airglow lines (particularly 
the Lyman-/? and Lyman-7 lines of HI) and using the relative strengths found in Feldman et 
al. (2001). These comparisons indicate that telluric NI emission is not the dominant cause 
of this zero-level problem. 

There are some cases where the 1134 A triplet is weak enough that the while the lines 
are intrinsically saturated, the observed profiles are not (i.e. the intensities at the center of 
the lines does not reach zero and the cores are not flat). In these cases the instrumental 
profile dominates the observed profile, and therefore Gaussian fits are justified and corrections 
cannot be made to the zero point of the spectrum. In the majority of cases, however, flat 
cores are observed and Voigt profiles are used. 

In the case of HD 53367, the line profiles are slightly assymetric and not well-fit by 
single- velocity component profiles, either Gaussian or Voigt. We measure the equivalent 
widths of the 1134 A triplet in this line of sight by a simple summation relative to the 
assumed continuum level, rather than the equivalent width of a fit. Errors are estimated 
based on the noise level in the continuum. 
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3.3. 1160 A Doublet 

The weak doublet at 1160 A has been the staple of NI studies utilizing Copernicus data. 
However, most HST spectra do not extend down to this wavelength. The one spectrum in 
our study that barely does (HD 147888) has inferior S/N in that region compared to the 
FUSE spectrum. Thus, this doublet is measured exclusively by FUSE. 

As shown in Table 4, the /-values for this doublet are significantly smaller than the rest 
of the absorption lines that we observe — the stronger component of the doublet has an /- 
value less than half of that of the next weakest line we observe (at 951.3 A), more than three 
orders of magnitude smaller than the weakest component of the 1134 A, and more than four 
orders of magnitude smaller than the weakest component of the 1200 A triplet. Therefore, 
when one or both components are detected, this doublet provides an important constraint 
on the overall curve of growth. Due to the small /-values of the doublet, Gaussian fits are 
always appropriate for both components (i.e. the observed line profiles are not saturated, 
show no resolved velocity structure, and the instrumental profile dominates even if there is 
some intrinsic saturation). In some spectra, a stellar feature interferes with the 1159.8 A 
component of the doublet, but in most cases the line can still be fit. 

There is some level of uncertainty in the /-values of these two lines (Morton 2003). The 
theoretically calculated line strengths cited in Tachiev & Froese Fischer (2002) differ from 
those in Hibbert et al. (1985) by approximately 0.068 dex (17%) for the stronger line and 
0.057 dex (14%) for the weaker line. Previous literature values show an even greater range; 
for instance, the theoretically calculated /-value of the 1159.8 A line in Cowan et al. (1982) is 
three times stronger than the astrophysically-derived empirical value in Lugger et al. (1978). 

For simplicity and self-consistency, we uniformly assume the Tachiev & Froese Fischer 
/-values. If our curves of growth are dominated by the weak lines, there will be a systematic 
difference between our derived column densities and those of Meyer et al. (1997) and Knauth 
et al. (2003). Since the Tachiev & Froese Fischer (2002) /-values for this weak doublet are 
larger than the Hibbert et al. (1985) /-values, we expect our derived column densities to be 
smaller if and when there are systematic differences. For simplicity, in this paper the values 
we quote from Meyer et al. and Knauth et al. are adjusted downward by 17%. 

3.4. 1200 A Triplet 

In the lines of sight with i^STdata, we are able to measure a triplet at 1200 A. The lines 
of this triplet are stronger than all of the other nitrogen lines measured in this study. Voigt 
profile fits are always used to fit these lines. There is a line of Mnll in the blue wing of the 
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bluest line of the triplet; this Mnll line is relatively weak and is removed with a simple fit (it 
is nearly Gaussian, convolved with our empirical STIS PSF). A stellar feature also interferes 
with this blue line of the triplet in some cases. However, we were able to reconstruct the 
local continuum and fit the line. 



3.5. Other Lines 

There are other neutral nitrogen lines that fall in the wavelength range covered by 
FUSE. However, none of these lines was clearly detected, due primarily to interference by 
strong H 2 absorption (see Table 5). There are no other UV absorption lines of nitrogen in 
the range covered by HST. 

3.6. Errors 

Errors on equivalent width measurements take two forms. When absorption lines are 
fitted to a Gaussian, errors are taken from standard error propagation of the curve-fitting 
routine and the functional form of a Gaussian. These errors include continuum placement 
error added in quadrature with the errors in the fit. When absorption lines are fitted to a 
Voigt profile, errors are not clear due to the lack of a functional form (the Voigt profile is fitted 
iteratively and involves convolution with the instrumental resolution element). Therefore, 
we estimate 1-a errors by first taking the difference between the fit and the data (both 
normalized) at each point. The uncertainty is then estimated to be the standard deviation of 
this difference, multiplied by the wavelength range of the fit. This method implicitly includes 
continuum placement errors. The validity of this method of error estimation was confirmed 
by using both methods on absorption lines for which a Gaussian was an appropriate fit. 

Potential errors from the placement of the zero level of the spectrum are not formally 
accounted for; however, the FUSE spectra of the lines of sight in this study have very wide 
H2 bands that constrain errors in the placement of the zero level over wide ranges of the 
FUSE bandwidth (Rachford et al. 2002); in most cases, the zero level is uncertain to only 
a few percent of the continuum level (with results in a correspondingly small uncertainty in 
equivalent widths due to this effect). We have commented above (§3.2) on the particular 
zero-level issues encountered above with the 1134 A triplet. 
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4. METHODS 

4.1. The Curve-of-Growth Method 

In order to derive nitrogen column densities, we follow the same procedure outlined in 
our previous work on 01 (Jensen et al. 2005). We fitted the measured equivalent widths 
of the various NI absorption lines to a curve of growth with a single velocity dispersion. 
Several curves of growth were constructed using the damping constants of the strongest 
lines (the lines of the 1200 A triplet, the 1134 A triplet, and the two observable lines of the 
953 A triplet), so that each absorption line was compared to the curve corresponding to its 
damping constant. We found the minimum \ 2 with respect to the two free parameters of 
column density N and 6-value, then used the extrema of the confidence ellipses to find the 
l-a errors in N and b- value. The best fits for these curves of growth are reported in Table 8, 
with l-a errors. Also reported in Table 8 are the abundances of nitrogen relative to hydrogen, 
given in parts per million (ppm). The errors on these nitrogen abundances are taken from 
error propagation done in standard (i.e. non- logarithmic) notation. In one case, HD 41117, 
the large errors in the NI column density combined with standard error propagation results 
in lower limit errors larger than the absolute value of the nitrogen abundance, and thus we 
have artificially capped the lower limit to be N/H > 1 ppm. 

There are potential systematic errors stemming from our assumption that the velocity 
structure of each line of sight is reasonably approximated by a single 6-value. However, 
when measurements from absorption lines with a wide range of /-values are available (our 
measurements typically span up to ~ 4 orders of magnitude in /-value), these systematic 
errors should be small. We discuss this idea in more detail in §4.4. 

Adopted curves of growth are shown in Figures 2-4. In the rest of this section, we 
describe the measurement of hydrogen column densities and the alternative methods and 
consistency checks that were employed in the final determination of nitrogen column densities 
and b- values for certain lines of sight. Again, we note here that these consistency checks rely 
on the assumption that the velocity structure is reasonably approximated by a single b- 
value. However, these consistency checks are only intended to rule out very disparate results 
(e.g. a column density from the curve of growth that implies a damped profile in one of 
the strong lines when such a profile is not observed) and not to differentiate between more 
subtle differences in column density and velocity structure. 
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4.2. Hydrogen Column Densities 

Molecular hydrogen column densities are taken from recent FUSE studies (Rachford et 
al. 2002; B. L. Rachford et al., in preparation). In these papers iV(H2) is determined for 
each sightline by fitting several low- J lines. Atomic hydrogen densities are taken from several 
sources, listed in Table 2. We adopt literature values where they exist. Table 2 summarizes 
this data. 

Bohlin et al. (1978) showed for less reddened lines of sight that there exists a correla- 
tion between selective reddening and total hydrogen column density such that iV(H tot ) = 
5.8 x l0 21 E B -v- Rachford et al. (2002) showed that this correlation remains valid to approx- 
imately ±0.30 dex (a factor of 2) in logA^(H to t) even for large hydrogen column densities 
(logiV(H tot ) > 21) and large molecular fractions of hydrogen; furthermore, the correlation 
is much more precise than ±0.30 dex in most cases. For two lines of sight (HD 38087 and 
HD 53367), we estimated total hydrogen column densities using this correlation, and then 
inferred atomic hydrogen column densities based on our molecular hydrogen results, that 
is, the atomic hydrogen column density is given by iV(HI) = 5.8 x 10 21 Eb~v ~ 2A^(H 2 ). 
We assume the error in iV(HI) using this method is ±0.30 dex, and the final error in the 
total hydrogen column density is determined from the error propagation in adding iV(H 2 ) 
and JV(HI); therefore lines of sight with a large molecular fraction and small error bars on 
iV(H 2 ) will have the smallest error bars in our total estimated iV(H tot ). Similarly, for HD 
179406 we quote the literature value of Hanson et al. (1992), which is not a direct mea- 
surement of HI but comes from the relationship found by Shull & Van Steenburg (1985), 
N(m) = 5.2 x 10 21 E B _ V . 

iV(HI) is determined in the remaining 27 lines of sight through profile fitting of Lyman- 
a. A few of these lines of sight warrant mentioning. First, we made our own measurement of 
the Lyman-a profile for the line of sight toward HD 40893, using STIS data sets O8NA02010 
and O8NA02020. Our analysis of the profile of the line yields a value for the column density 
of log A r (HI) = 21.50 ± 0.10. This result is somewhat larger (by 0.22 dex) than the result 
that would be obtained through use of the Bohlin et al. (1978) relationship, but within our 
estimated limits of that relationship. 

Second, in lines of sight studied by Diplas & Savage (1994), stellar contributions to the 
Lyman-a profile were estimated. Of the stars where we cite interstellar HI along the line of 
sight from Diplas & Savage, only HD 42087 is B2.5 or cooler. Diplas & Savage deemed the 
estimated stellar contamination to be insignificant. 

Third, we quote the values of Cartledge et al. (2004) for HD 27778 and HD 147888. 
These stars are cool enough that stellar Lyman-a might be a concern. For HD 27778, 
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however, the derived iV(HI) matches, within the errors, the value estimated by Rachford et 
al. (2002) using the Bohlin et al. (1978) relationship. It also results in an interstellar Kr/H 
abundance that precisely matches the interstellar mean. For HD 147888, Cartledge et al. 
used HD 147933 as a proxy. However, they report that the Lyman-a profile of HD 147888 
results in the same column density. This column density also agrees, within the errors, with 
an earlier literature value of Fitzpatrick & Massa (1990). Thus, for HD 27778 and HD 147888 
we conclude that stellar contamination is not significant. 



4.3. Alternate Methods 

For each sightline, we use column density and 6-value results from the curve-of-growth 
analysis to generate profiles of multiple absorption lines (in particular the 1134.98 A line and 
the 1160 A doublet) as a check for consistency. These profile fits occasionally reveal flaws 
in the curve-of-growth analysis. In these cases, we use other methods to determine column 
density. One method we used was to assign a uniform fractional error (based on the average 
of the fractional errors) to each of the equivalent width measurements, then analyze the curve 
of growth. This method tends to put more emphasis on the weak lines and less emphasis on 
the strong lines, which may be justified if there is systematic error in the strong lines that 
is unaccounted for. We use this method for three sightlines (HD 24534, HD 42087, and HD 
185418). In each of these three cases, the "uniform error" curve-of-growth method produces a 
result which is more consistent with line profiles than the standard curve-of-growth method. 

Occasionally, the x 2 array from our curve-of-growth analysis has multiple local minima. 
In these cases, we must select which solution is most consistent with the observed profiles 
of the lines. This is the case for HD 24534 and HD 53367. In both cases we select as our 
solution the point which is a local minimum, but not the global minimum. In the case of HD 
24534, this is done for the "uniform error" curve-of-growth result. The selected solutions 
and quoted errors are allowed to within 1-a, while the other local minima are excluded on 
the basis of inconsistency with the observed profiles of the absorption lines. 

The fact that curve-of-growth method produces two local minima indicates complex 
velocity structure. For HD 24534, we noted the hints of complex velocity structure for OI 
(Jensen et al. 2005); this is even more clear in high- resolution optical KI and Cal data 
(Welty k Hobbs 2001; Welty et al. 2003; Pan et al. 2004). We also noted above (§3.2) that 
the profiles of lines in the 1134 A triplet were assymetric in the case of HD 53367. 
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4.4. The Linear Approximation and Saturation Concerns 

For each line of sight we calculated the resulting column density if the weak lines of the 
1160 A doublet are assumed to be on the linear portion of the curve of growth. In one line 
of sight (HD 186994) we are unable to measure either component of the 1160 A doublet. 
Of the remaining 29 lines of sight, 28 show consistency within the 1-a errors between the 
column density results of the curve-of-growth method and the weak-line assumption for the 
weakest line measured (the 1160.9 A line when measured, otherwise the 1159.8 A line). The 
other line of sight, HD 197512, is consistent to 2-cr. This demonstrates consistency in that it 
shows that the strong lines do not dominate the determination of column densities, and that 
results from using our curve-of-growth method would be similar to the results of a study if 
we used the weak lines only, assuming them to be unsaturated. 

Systematic errors might still be a problem if the velocity structure cannot be reasonably 
approximated by a single-component velocity dispersion. However, the internal consistency 
of most of our curves of growth is what leads us to believe that this is not the case. When 
both components of the weak doublet are detected, the ratio of equivalent widths is almost 
always very nearly what would be from expected from assuming both lines are unsaturated 
(~ 3.6), implying little to no saturation. For the cases where only the stronger component of 
the doublet is expected, the lack of a detection of the weaker line indicates that the stronger 
line cannot be significantly saturated. The case where the greatest saturation (as determined 
from the 1160 A doublet ratio) is detected is for HD 24534, which as noted above, has a 
very complex velocity structure. However, our final derived column density still places both 
components of the 1160 A doublet near the linear portion of the curve, and our result is 
consistent (within the errors) with the column density derived by Knauth et al. (2003) based 
on the weak doublet only. 

To the extent that our measured sets of equivalent widths for some lines of sight do 
deviate somewhat from a single velocity-dispersion curve of growth (e.g. note the "hitch" in 
the data points near the linear/flat transition for the curve of growth of HD 210839), we can 
still be confident that the correct column density has been measured if the 1160 A doublet 
ratio is similar to the expectation from the unsaturated assumption. This is true for all 
cases where both components of the doublet are measured except HD 24534. In any case, 
determining complex velocity structure solely from FUSE data would be very imprecise at 
best, and we cannot do a more complex (i.e. multiple-component) curve-of-growth analysis 
without a priori assumptions about the velocity structure. High-resolution optical data on KI 
and Cal is available for some of the lines of sight in our sample (Welty & Hobbs 2001; Welty 
et al. 2003; Pan et al. 2004), but assuming that the relative abundances of these elements to 
NI are constant from cloud to cloud in a given line of sight or that b- values are the same for 
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all elements in a given cloud (especially given the different ionization potentials) is not more 
justified than assuming that the entire velocity structure is reasonably approximated by a 
single velocity dispersion and/or that the weak 1160 A doublet is not significantly saturated. 

5. RESULTS AND DISCUSSION 

Adopted curves of growth, with column density and b- value results, are shown in Figures 
2-4. We find an error- weighted average for the N/H ratio in the ISM of 51 ± 4 ppm (error in 
the mean). This is smaller than the average N/H abundance found by Meyer et al. (1997) of 
62^3 ppm. The discrepancy is more than 1-a but within 2-a. On the other hand, our results 
are in nearly perfect agreement to the weighted average of Knauth et al. (2003), 51 ±3 ppm. 
There are a handful of lines of sight that overlap between our study and Knauth et al. (see 
§5.2); the weighted average of the lines of sight that are unique to the Knauth et al. study 
is 53 ± 4. Again, both of these weighted averages have been adjusted downward by 17% to 
account for the systematic difference in /-values. 

A plot of N/H ratios as a function of total hydrogen column density is shown in Figure 
5 — lines of sight from Knauth et al., Meyer et al., Jenkins et al. (1999), and Sonneborn et 
al. (2000) are included for comparison, although we do not include the Jenkins et al. and 
Sonneborn et al. results in our correlation or weighted average calculations. Attempting 
both linear regressions and a four-parameter Boltzmann function similar to that used by 
Cartledge et al. (2004, 2006), we do not find statistically significant correlations of the N/H 
ratio with regard to various extinction parameters (total extinction Ay, selective extinction 
Eb-Vi or the total-to-selective ratio Ry) or total hydrogen volume density (nn) — all of these 
calculations produce fits that are consistent with a flat line (i.e. no correlation) to within 
approximately 1-a. 

Two statistically significant correlations do come to light. First, a simple linear regres- 
sion of N/H vs. N(H tot ) in our data shows a negative slope that is significant to approxi- 
mately 3-a. When we add in the lines of sight from Meyer et al. and the unique lines of 
sight from Knauth et al. (i.e. the ones that we have not independently analyzed), this anti- 
correlation remains, with approximately the same value of the slope and the same statistical 
significance. However, by excluding HD 152236 (a line of sight with an unusual value and a 
small error; see §5.1), the significance of the correlation is significantly lessened. 

We also analyzed this potential correlation with a four-parameter Boltzmann function of 
the form y = A H — B ^ A XQ , where in this case y is the nitrogen abundance, x is the hydrogen 

l+e m 

column density, x is a transition value in the hydrogen column density, m is the local slope 
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of the transition, and A and B are the asymptotic values of y for high and low column 
density, respectively. Using this formula, the best fit, with a reduced x 2 °f 0.138, has an 
abrupt transition at log JV(H tot ) = 21.05 with values of 48 ppm and 60 ppm for the high and 
low column density limits, respectively. Despite the small x 2 , however, the errors in the two 
asymptotic values overlap. 

The second correlation is with the molecular fraction of hydrogen, /(H 2 ). We observe a 
positive correlation between N/H and /(H 2 ) with a slope statistically significant to approx- 
imately 3-0" based merely on our own data. However, this correlation disappears when we 
add in the Meyer et al. and Knauth et al. results, and can mostly likely be disregarded. 

In an attempt to further investigate the potential correlation with total hydrogen column 
density, we analyzed previous Copernicus equivalent width measurements from Bohlin et 
al. (1983), with the goal of exploring lower column density lines of sight (these are the same 
equivalent widths analyzed by York et al. 1983). To avoid systematic errors, we used the 
same /-values and curve-of-growth methods that we used to analyze our sample. However, 
the results showed a significant amount of scatter, and there are some potential systematic 
errors, such as interference from stellar lines. Additionally, some of the equivalent width 
measurements in certain lines of sight in the Bohlin et al. (1983) data are not self-consistent, 
with some weak lines measured to have greater equivalent widths than other lines with 
stronger /-values. The table in Bohlin et al. that probes smaller column densities is Table 4. 
However, Table 4 only utilizes a triplet of NI at 952 A. Our analysis of the lines of sight in this 
table do show small abundances relative to hydrogen (a weighted average of ~ 20 ppm for 10 
lines of sight). However, there are potential systematic errors for two reasons. First, there 
may be errors in the equivalent width measurements, as Copernicus suffered from significant 
scattered light in the region below 1000 A. Scattered light would systematically reduce the 
equivalent widths, so it is not surprising that these results are small. Secondly, because we 
did not use the 952 A triplet in our more reddened lines of sight (due to interference from 
H2), we do not have any examples in our data of the self-consistency (or lack thereof) of 
the 952 A triplet with the lines we did use. Any systematic error between the 952 A triplet 
/-values and the /-values of the lines used in this paper would lead to similar systematic 
error between the final results for the two samples. 

Many other sources in the literature are also not easily compared to our results for a two 
reasons. First, individual equivalent widths, column densities, and/or abundances are not 
always explicitly given. Secondly, as with Table 4 of Bohlin et al. (1983), different absorption 
lines are measured due to the differences in column density. One example of this is Lugger 
et al. (1978), who measured the nitrogen content in 10 lines of sight. Though Lugger et al. 
(1978) used the 1134 A and 1200 A triplets, they did not measure the same weak lines as 



-16- 



we did in this study except in a few cases. Even the additional literature values (Jenkins et 
al. 1999; Sonneborn et al. 2000) cited in Knauth et al. do not meet this criterion, having 
measured lines such as the 952 A triplet (as well as the 1134 A triplet, though mostly to 
analyze structue). The Meyer et al. sample, although somewhat small, is the best sample 
that meets the following criteria: (1) smaller column densities (log-/V(H to t) ~ 20 — 21) that 
are not too small such that ionization might be significant and (2) measurement of at least 
some of the same weak lines. Lines of sight with even smaller column densities are generally 
shorter lines of sight, where the effects of ionization in the LISM/Local Bubble become 
significant. 

Knauth et al. concluded that nitrogen was depleted in lines of sight with total hydro- 
gen column density of logiV(H tot ) > 21. This was based on a student's t-test; Knauth et 
al. calculated the probability that the two samples (the logA^(H to t) < 21 sample and the 
logiV(Htot) > 21 sample) come from the same parent sample was 7.6%. This probability 
was lessened to 1.3% by the exclusion of the outlying point HD 36486 (5 Ori), from Jenkins 
et al. (1999). However, the standard form of the t-test compares the standard deviations of 
the two samples but does not take into account errors in the individual points. Knauth et al. 
do not explicitly state if and how they take the individual errors into account. Similar t-test 
calculations with our sample fail to yield results similar to those of Knauth et al.. However, 
it is entirely possible this is an artifact of the broader range of errors found in our sample. 
Since the t-test does not account for these errors, lines of sight in our sample that deviate 
from the mean but also have a large error are weighted more in the t-test than is justified. 

Given that the t-test is not appropriate for our sample, in addition to the correlation fits 
discussed above, we examined the weighted averages of N/H in different regimes of hydrogen 
column density. We have already discussed above the differences in our weighted average of 
N/H vs. that of Meyer et al. and Knauth et al. — the weighted averages for our results and 
those of Knauth et al. differ from the weighted average of the Meyer et al. lines of sight by 
more than l-cr. The difference in iV(H to t) between these samples is significant; 29 of our 30 
lines of sight have more hydrogen than the highest column density lines of sight in Meyer et 
al. (HD 143275 and HD 149757); 22 of our 30 lines of sight have a hydrogen column density 
that is at least a factor of two greater. Combining the samples of Meyer et al., Knauth et 
al., and this paper (using our results for lines of sight common to this paper and Knauth 
et al.), we compare the average N/H in lines of sight of line with log-/V(H to t) < 21 against 
those with logA^(H to t) > 21; the weighted averages overlap slightly within the errors. If we 
increase the breaking point to logiV(H tot ) = 21.1, the averages cease to overlap within the 
errors. At a cutoff of logiV(H tot ) = 21.4, which roughly halves the combined sample, the 
weighted average is 59 ±3 ppm for lines of sight with log iV(H t ot) < 21.4, and 49 ±4 ppm for 
lines of sight with logA^(H tot ) > 21.4. 
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We have not included the Jenkins et al. (1999) and Sonneborn et al. (2000) results in 
these fits and weighted average calculations because they are measured using largely different 
absorption lines, as discussed above. If we do include those results, these correlations and 
weighted averages generally remain unaffected; if we include the two Sonneborn et al. (2000) 
lines of sight but exclude the outlying point HD 36486 (as Knauth et al. did to show increased 
statistical significance), the correlations and differences between the weighted averages are 
strengthened somewhat. 

Based on these weighted averages, combined with the anticorrelation with iV(H tot ) dis- 
cussed above, we conclude that that we are seeing some evidence, however limited, that 
nitrogen depletion does increase in lines of sight with total hydrogen larger column densities. 
This supports the conclusion of Knauth et al., but in both cases there are reasons for caution. 
Our statistical measures are greater than 1-a but in some cases not by much. Furthermore, 
one line of sight that influences these results is HD 152236. As discussed in §5.1, this line of 
sight may have an unusually low N/H ratio due to nucleosynthetic effects rather than cloud 
conditions and chemistry (incorporation into grain material, direct depletion onto grains, 
and/or incorporation into nitrogen-bearing molecules). We will proceed in this paper with 
the tentative conclusion that nitrogen is mildly depleted in lines of sight with greater column 
density — noting the uncertainties — and discuss possible reasons for this variation, as well as 
some side issues. 



5.1. HD 152236 

Most of the individual results we adopt are not substantially different from the previously 
measured abundance of nitrogen in the interstellar medium. The primary exception is our 
result for HD 152236. Our curve-of-growth analysis measures an N/H ratio of 26^g ppm 
(several a below the Meyer et al. mean, even after accounting for any systematic differences 
in /-values) and a 6-value of 15.9 ± 1.4 km s -1 . While this is an unusually large 6-value, 
it is consistent with what we previously found for OI in this line of sight (Jensen et al. 
2005), as well as the line profiles of both NI and OI. What is unusual, however, is that the 
measured column density for OI in our previous paper was very large, leading to the largest 
overall O/H ratio in that sample, while here the measured N/H ratio is the second smallest 
in our sample (the line of sight toward HD 37903 has a smaller N/H ratio but with a much 
larger upper limit). The consistency of the b- values leads us to believe that even if there is 
systematic error in the N/H or O/H ratios, there is unlikely to be systematic error in NI 
column density relative to OI or vice versa, and the N/O ratio must be very small. Using 
the adopted numbers in this paper and Jensen et al. (2005), N/O = 1.2 x 10~ 2 , compared 
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to the solar (Holweger 2001) and ISM (Knauth et al. 2003) ratios where N/O - 0.2. Even if 
a smaller column density for OI is accepted (as discussed in Jensen et al.), N/O is unlikely 
to be any larger than 0.05. 

It is unclear what process or processes in the ISM could cause such a low N/O ratio. 
To explain stellar abundances, Walborn (1976) suggested that some members of the Sco 
OBI association, of which HD 152236 is a member, may have formed from nitrogen-deficient 
material; this result may be evidence for that suggestion. Despite this piece of corroborating 
evidence, the ISM processes or chemistry that may have initially created this enhanced 
depletion are not immediately obvious. Furthermore, it is not known to what extent the 
abundances of an OB star trace the abundances of the intervening ISM. Whether HD 152336, 
with the largest total hydrogen column density in our sample, shows such a low N/H ratio 
due to nucleosynthetic effects or cloud conditions and chemistry is ultimately unclear. 

5.2. Comparison With Previous Work 

Our sample overlaps with the sample of Knauth et al. (2003) for eight lines of sight: HD 
24534, HD 73882, HD 110432, HD 147888, HD 179406, HD 185418, HD 192639, HD 210839. 
Our final adopted N/H ratios agree with those of Knauth et al. for all eight sightlines within 
the errors, even after potential systematic errors from differing /-values are considered. It is 
worth noting that our method of adopting of column densities for HD 24534 and HD 185418 
was different than for the other sightlines in this sample (see §4.3). In general, our error 
bars are larger than those found by Knauth et al.; this is likely an artifact of our use of the 
curve-of-growth method as opposed to the determination of column densities from the weak 
doublet only. 

5.3. b- values 

It has long been assumed (York et al. 1983) that interstellar OI and NI are coincident 
in position and velocity space. If this is true, the two elements should have the same b- value, 
and if a 6-value can be measured for one element that same 6-value can be assumed for the 
other element. Using the same method and data sets as our previous work on OI (Jensen et 
al. 2005) allows us to compare our results for the independently measured b- values of OI and 
NI. Of the 24 common sightlines for which we report both OI and NI column densities (we 
did not report OI results for HD 40893, HD 41117, HD 147888, HD 149404, HD 167791, and 
HD 203938), we find 1-a agreement for b- values for 15 lines of sight. Another six lines of sight 
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(HD 24534, HD 37903, HD 38087, HD 154368, HD 170740, and HD 206267) fall outside the 
region of 1-a agreement but agree within 2-<r. Three lines of sight (HD 38087, HD 199579, 
and HD 207198) only possess agreement within 3-a. Therefore, with the possible exception 
of a few cases, disagreement in b- values is likely statistical. 

The three lines of sight with only 3-a agreement are the most likely to be cases of either 
systematic errors somewhere in the measurement of absorption lines (affecting the curve of 
growth and subsequent calculation of the 6-value) or real cases where 01 and NI are not in 
the same clouds and/or do not have the same physical conditions. In the case of HD 38087, 
the column density and 6-value for OI we derived in Jensen et al. (2005) are very uncertain, 
based on only two absorption lines. It is probable that these values are more uncertain than 
the quoted errors. Therefore, the discrepancy in 6-values for NI and OI for HD 38087 cannot 
be considered significant. 

HD 199579 was flagged as a potential outlier in Jensen et al. (2005). A very large 
b- value and a very small oxygen column density (leading to a small O/H ratio) were mea- 
sured. However, profile fitting of the 1039 A absorption line of OI seemed to favor a larger 
column density and smaller 6-value of ~ 10 km s" 1 — consistent with our 6-value measured 
for nitrogen. Since our measurement of NI for HD 199579 does favor this 6-value and does 
have weak lines constraining the linear portion of the curve of growth, it is probably best to 
adopt the NI b- value. 

The situation is somewhat more complex for HD 207198. The NI curve of growth is 
relatively straightforward, with several strong lines constraining the flat portion of the curve, 
and several weak lines constraining the linear portion and the transition region between linear 
and flat. The OI curve of growth is somewhat less certain visually, though the statistical 
analysis constricts the range of possible 6-values. In Pan et al. (2004), optical data on KI 
shows multiple cloud components. However, the cloud components are roughly symmetric 
and heavily blended. The total KI profile (i.e. approximating the blended profiles of multiple 
cloud components as one profile) has a FWHM of ~ 12kms _1 , which would correspond to 
a b- value of ~ 7kms~ 1 . If a b- value is determined empirically for each absorption line of 
OI in the HD 207198 spectrum, these empirical 6-values increase with line strength. This 
indicates a complex line of sight with small clouds that become important in the stronger 
lines. It would be difficult to perform meaningful analysis of this line of sight using a multiple- 
component profile, due to the limits of the resolution of FUSE (~ 15 — 20 km s _1 ). However, 
the 6-value of all cloud components taken as a whole is probably better represented by the 
~ 7 km s _1 6-value that we estimate from Pan et al. (2004) and our NI results than our 
previous results for OI (Jensen et al. 2005). 
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5.4. N/H of Galactic Stars 

The two current solar N/H ratios that are widely accepted are 93 ± 16ppm, from 
Grevesse & Noels (1993), and the more recent work of Holweger (2001), who found (N/H) = 
85 ± 22 ppm. As with many other elements in the ISM, however, the solar value may not 
be an appropriate abundance standard. Lodders (2003), for example, suggests a pre-solar 
nebula N/H abundance of 79+Jg ppm. N/H abundances found in B stars further complicate 
the situation: Snow & Witt (1996) find average stellar N/H ratios of 63 ppm (field B stars) 
and 68 ppm (cluster B stars), while Sofia & Meyer (2001) find an average N/H of 65 ppm. 
As discussed in Meyer et al. (1997), it is therefore likely that B star abundances are a better 
standard for N/H in the ISM than the solar abundance. These surveys did not measure N/H 
in the same samples of F and G stars where other abundances were measured. Therefore, 
the N/H abundance of B stars would seem represent a better potential ISM standard than 
the solar N/H abundance; however, we do not have the full picture, as it is unclear how the 
N/H abundance of F and G stars measures up. It is also worth noting that B stars seem to 
be a worse potential abundance standard than F and G stars for many other elements (Sofia 
& Meyer 2001). 

If the N/H abundances of B stars are adopted as the ISM standard, the situation is 
rectified somewhat: our results would indicate that 14 ppm of nitrogen is missing from the 
gas-phase ISM, without taking into account the errors. The results of Knauth et al. (2003) 
would indicate a similar gas-phase nitrogen deficiency. On the other hand, the lines of sight 
in Meyer et al. appear to have very little gas-phase NI depletion. However, the errors 
in B star abundances as determined by Sofia & Meyer (2001) are quite large (34 ppm). 
Additionally, these inferred depletions change with differing /-values for the weak lines, as 
we have noted many times already (the derived stellar N/H abundances are not dependent 
on these /-values) . It is therefore difficult to draw any firm conclusions regarding the amount 
of nitrogen in dust, even if a B star abundance is adopted. Even if the absolute amount 
is not known, however, the difference in N/H of ~ 10 ppm for lines of sight with larger 
hydrogen column densities that our results hint at must be explained. We briefly discuss a 
few possibilities in §5.6. 

As discussed in many previous papers (e.g., Snow 2000; Sofia & Meyer 2001; Jensen et al. 
2005), it may be the case that good abundance standards in the ISM do not exist. Processes 
that occur during the formation of stars such as ambipolar diffusion, photoevaporation of 
the disk, and late metallicity pollution from the infall of protoplanets are likely to have a 
great impact on the details of each individual star, explaining why the variation in stellar 
abundances is greater than the variation of the well-mixed ISM from which they formed. If 
a better understanding of an appropriate ISM standard for N/H can be gained, however, we 
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may be able to better understand just where the NI goes in the densest sightlines. 

5.5. Ionization Effects 

It is fair to ask whether or not any of our results, or the results of previous studies, may 
be influenced by ionization effects. Indeed, Moos et al. (2002) concluded that ionization is 
responsible for a deficiency of NI relative to HI in the local interstellar. However, Jenkins et 
al. (2000) discuss the processes that couple the ionization fractions of nitrogen and oxygen 
to hydrogen. With a larger ionization cross-section, it might be expected that nitrogen is 
more ionized than hydrogen. However, the ionization fractions of nitrogen and hydrogen are 
also governed by a charge-exchange reaction (Butler & Dalgarno 1979). Only if n c ^ nn will 
nitrogen and hydrogen show a significant departure from the very similar ionizations assumed 
in §1. Figure 2 of Jenkins et al. (2000) shows a model that indicates that nitrogen will be 
substantially more ionized than hydrogen only if log iV(HI) <C 19. However, in lines of sight 
toward four white dwarfs Jenkins et al. did find a smaller abundance of NI than expected. 
This indicates potential problems with the model, possibly because high-energy photons 
from helium recombination are a large influence (the model assumes only photoionization 
from hot stars and a hot-gas conduction front). 

If there is a potential problem with the Jenkins et al. (2000) model, it is not clear how 
this should affect the assumption that ionization is not important in our lines of sight. The 
lines of sight in this study have hydrogen column densities at least two orders of magnitude 
larger than the four lines of sight in the Jenkins et al. study; the column densities in the 
Meyer et al. (1997) lines of sight are at least one order magnitude larger. The reddening of 
our lines of sight also imply fairly large dust content, which might contribute to shielding the 
gas clouds from higher-energy photons. Our lines of sight are also not necessarily subject to 
the unique conditions of the LISM/Local Bubble. 

An additional consideration is that all of the lines of sight in this study have /(H 2 ) > 
0.09, including nearly half with /(H 2 ) > 0.4. Lines of sight with such large molecular 
fractions of hydrogen are typically indicative of environments that are largely self-shielded. 
In seven lines of sight spanning several orders of magnitude in /(H 2 ), Meyer et al. found no 
evidence of increasing depletion of nitrogen as the molecular fraction of hydrogen increases 
(see Figure 2 of Meyer et al. 1997). In fact, the same figure shows very limited evidence 
that nitrogen is more depleted in lines of sight with smaller /(H 2 ) — this would be consistent 
with the possibility that ionization effects are important. N/H follows a similar pattern with 
respect to total hydrogen column density as well — the two lines of sight with the smallest 
N/H also have two of the three smallest hydrogen column densities. However, it is worth 
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noting that the Meyer et al. sample, at a total of seven lines of sight, is very small. 

5.6. Dust Grains and Sinks of NI 

Given that the results of this paper and Knauth et al. (2003) point (however tentatively) 
toward an increase in nitrogen depletion in dense lines of sight, the question must be asked: 
where does the NI go? Knauth et al. determined that N 2 could not be a significant sink for NI 
in the lines of sight in that study. Furthermore, unlike oxygen, silicon, and iron, nitrogen is 
not a significant ingredient in most dust models — a moderate amount of nitrogen is required 
for models that include organic refractory material (Snow & Witt 1996), but nitrogen is 
not necessary for dust models dominated by silicates and oxides. Perhaps, as Knauth et 
al. suggest, our understanding of interstellar nitrogen chemistry is incomplete. Another 
possibility is that our observations are probing lines of sight with one or more regions where 
the local density is high enough to see a fair amount of atomic nitrogen depleted directly 
onto dust grains, as opposed to being an ingredient in the molecules that make up grains. 
In this case no further understanding of nitrogen chemistry is needed. 

6. SUMMARY 

We measured the NI column densities and b- values of 30 lines of sight with log A(H tot ) > 
21 using measurements of several NI absorption lines and fitting those to curves of growth. 
We find a weighted-average value for N/H in these sightlines of 51 ± 4 ppm. Within the 
uncertainties, this is similar to the /-value adjusted ratio of 51 ± 3 found by Knauth et al. 
(2003), but smaller than the /-value adjusted ratio of 62^ found by Meyer et al. using less 
reddened lines of sight. We also note a trend with respect to A(H tot ) when all three samples 
are considered; similar statistically significant trends do not exist for extinction parameters 
or other parameters such as /(H2) and nn- We take these results to be limited evidence that 
the conclusions of Knauth et al. are correct, that the gas-phase depletion of nitrogen relative 
to hydrogen is enhanced in lines of sight where log A(H to t) > 21. Potential sinks of NI in 
the ISM are unclear. The reason for nitrogen depletion may be the direct accumulation of 
nitrogen onto large grains — an indication that perhaps the lines of sight observed are nearing 
the physical characteristics of the ever-elusive "translucent" clouds. Continued improvement 
in observational capability (such as with the Cosmic Origins Spectrograph) may allow these 
clouds to finally be detected. Finally, instead of explanations that invoke dust models or 
nitrogen chemistry, we cannot rule out nucleosynthetic effects as playing a significant role in 
these variations. Our prime example of this possibility is HD 152236, where low N/H and 
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N/0 ratios in the line of sight may be correlated (as either a cause or an effect) with similar 
ratios in the Sco OBI association. 
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Sample Spectra for HD 210839 




Wavelength (A) 



Fig. I. — The sightline toward HD 210839 has relatively high S/N, and all absorption lines 
used in this study are detected. Spectra are shown unnormalized and without velocity 
correction. Top Left: — The doublet at 951 A is shown, from SiC2A FUSE data. The two 
lines of the doublet are between 951.0 A and 951.5 A. The strong feature below 951.0 A is a 
blend of an H 2 line and an 01 line. Top Right: — The two visible components of a triplet at 
953 A are shown, from SiC2A FUSE data. The reddest component is destroyed by a strong 
H 2 feature. Middle Left: — The visible component of a doublet at 959 A is shown, from 
SiC2A FUSE data. The line is at approximately 959.5 A. The weaker, redder component 
may be seen at approximately 960.2 A but its separation from an H 2 line is unclear. The 
other features are H 2 lines. Middle Right: — The strong triplet at 1134 A is shown, from 
LiFIB FUSE data. The weaker feature at 1133.6 A is an Fell line. Bottom Left: — The 
weak doublet at 1160 A is shown, from LiFIB FUSE data. Bottom Right: — The strong 
triplet at 1200 A is shown, from HST STIS data. An Mnll line resides in the blue wing of 
the bluest component of the triplet. 
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Curves of Growth, HP 24534 through HP 110432 
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Fig. 2. — Curves of growth for HD 24534 through HD 110432. Measured equivalent widths 
of NI absorption lines are plotted on top of appropriate curves of growth for the adopted 
solution of column density and b- value. Equivalent width measurements are red squares with 
blue error bars. The solid, dotted, and dotted-dashed lines in the damping portion of the 
curve represent damping constants of 1.51 x 10 8 , 2.19 x 10 8 , and 4.07 x 10 8 , respectively. 
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Curves of Growth, HP 147888 through HP 197512 
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Fig. 3. — Curves of growth for HD 147888 through HD 197512. Measured equivalent widths 
of NI absorption lines are plotted on top of appropriate curves of growth for the adopted 
solution of column density and b- value. Equivalent width measurements are red squares with 
blue error bars. The solid, dotted, and dotted-dashed lines in the damping portion of the 
curve represent damping constants of 1.51 x 10 8 , 2.19 x 10 8 , and 4.07 x 10 8 , respectively. 
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Fig. 4. — Curves of growth for HD 199579 through HD 210839. Measured equivalent widths 
of NI absorption lines are plotted on top of appropriate curves of growth for the adopted 
solution of column density and 6-value. Equivalent width measurements are red squares with 
blue error bars. The solid, dotted, and dotted-dashed lines in the damping portion of the 
curve represent damping constants of 1.51 x 10 8 , 2.19 x 10 8 , and 4.07 x 10 8 , respectively. 
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og (N I / H) vs. Total Hydrogen Column Density 




-5.0 , , , , i , , I i , , , LU , , , I 

19.5 20.0 20.5 21.0 21.5 22.0 

log /V(Htot) WHtot) in cm" 2 ) 

Fig. 5. — A plot of log (N/H) as a function of log AT(H tot ). Blue Triangles: — Lines of 
sight from Meyer et al. (1997). Red Squares: — Lines of sight from Knauth et al. (2003) 
not analyzed in this study. Green Circles: — Lines of sight from the literature (Jenkins et 
al. 1999; Sonneborn et al. 2000) measured by IMAPS, included in the analysis of Knauth 
et al.. Black Diamonds: — Lines of sight from this study. Dashed Line: — The mean 
interstellar N/H found by Meyer et al., 62 ppm, is represented by the dotted line. This is 
very similar to the average N/H abundance of B stars found by Snow & Witt (1996) and 
Sofia & Meyer (2001) discussed in §5.4. Dashed-Dotted Line: — The average interstellar 
N/H of 51 ppm found by Knauth et al. and this study for lines of sight with larger hydrogen 
column densities. The point with the largest N/H ratio is HD 41117, which is not very well 
determined. The plotted point essentially represents an upper limit. Incorporating points 
from all three studies (but discarding the IMAPS points), we find a linear anticorrelation 
between N/H and iV(H tot ) of the form y = (—0.36 ± 0.13)x + (63 ± 3) where y is N/H in 
parts per million (ppm) and x is iV(H tot )/10 20 . See §5 and its accompanying subsections for 
further discussion of these results. 
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Table 1. Sightlines: Stellar Data 



Star Name Distance (pc) Distance Reference 11 1 


b 


Spectral Type 


HD 24534 


590 


2 


163.08 


-17.14 


09.5pe 


HD 27778 


223 


1 


172.76 


-17.39 


B3V 


HD 37903 


910 


2 


206.85 


-16.54 


B1.5V 


HD 38087 


480 


2 


207.07 


-16.26 


B5V 


HD 40893 


2800 


2 


180.09 


+4.34 


BOIV 


HD 41117 


1000 


2 


189.69 


-0.86 


B2Iae 


HD 42087 


1200 


2 


187.75 


+1.77 


B2.5Ibe 


HD 46056 


2300 


2 


206.34 


-2.25 


08V 


HD 46202 


2000 


2 


203.31 


-2.00 


09V 


HD 53367 


780 


2 


223.71 


-1.90 


BOIVe 


HD 73882 


1100 


2 


260.18 


+0.64 


08V 


HD 110432 


301 


1 


301.96 


-0.20 


B2pe 


HD 147888 


136 


1 


353.65 


17.71 


B5V 


HD 149404 


820 


2 


340.54 


+3.01 


09Iae 


HD 152236 


1800 


3 


343.03 


+0.87 


BHa+pc 


HD 154368 


960 


2 


349.97 


+3.22 


091a 


HD 167971 


730 


2 


18.25 


+1.68 


08e 


HD 168076 


2100 


2 


16.94 


+0.84 


05f 


HD 170740 


213 


1 


21.06 


-0.53 


B2V 


HD 179406 


160 


2 


28.23 


-8.31 


B3V 


HD 185418 


950 


2 


53.6 


-2.17 


B0.5V 


HD 186994 


2500 


2 


78.62 


+ 10.06 


BOIII 


HD 192639 


1100 


2 


74.90 


+1.48 


08e 


HD 197512 


1700 


2 


87.89 


+4.63 


B1V 


HD 199579 


1200 


2 


87.50 


-0.30 


06Vc 


HD 203938 


700 


2 


90.56 


-2.23 


B0.5IV 


HD 206267 


1000 


3 


99.29 


+3.74 


06f 


HD 207198 


1000 


2 


103.14 


+6.99 


09IIIe 


HD 207538 


880 


2 


101.60 


+4.67 


BOV 


HD 210839 


505 


1 


103.83 


+2.64 


06If 



a References: — (1) Hipparcos parallax of 4-ct precision or better. (2) Spectroscopic 
distance modulus. (3) Member of an OB association, cluster, or multiple-star system. 
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Table 2. Sightlines: Hydrogen Data a 



Sightlinc logJV(H 2 ) Rcf. b logiV(HI) HI Method c Ref. d logJV(H to t) n H /(H 2 ) 



HD 24534 


20 


92 


± 


0. 


.04 


1 


20 


.73 


± 





.06 


Ly-a 


1 


21.34 ± 0.03 


1 


.21 





76 


HD 27778° 


20. 


79 


± 


0. 


06 


1 


21. 


10 


± 





10 


Ly-a 


2 


21.40 ± 0.06 


3 


.62 





.49 


HD 37903 


20. 


92 


± 


0. 


06 


2 


21 


.17 


± 





.10 


Ly-a 


1 


n-i ™ + 0.06 

21 -50Io.o5 


1 


.12 





.53 


HD 38087 


20. 


64 


± 


0. 


07 


2 


20. 


.91 


± 





.30 


Bohlin 


3 


21 - 23 lo.l3 


1 


.14 





.52 


HD 40893 


20. 


.58 


± 


0. 


05 


2 


21 


.50 


± 





.10 


Ly-a 


3 


21.59+^.og 





.49 





.19 


HD 41117 


20. 


68 


± 


0. 


10 


2 


21. 


40 


± 





.15 


Ly-a 


1 


o-i ^,,+0.12 

21.54^ n 


1 


.12 





,28 


HD 42087 


20. 


.52 


± 





12 


2 


21 


.39 


± 





11 


Ly-a 


1 


21.49 ± 0.09 





.84 





.21 


TTT~\ A t~T\T" 

HD 46056 


20. 


68 


± 


0. 


.06 


2 


21 


38 


± 





.14 


Ly-a 


1 


91 to + 0.11 





.47 





,29 


HD 46202 


20. 


68 


± 


0. 


07 


2 


21 


.58 


± 





15 


Ly-a 


1 


21.68 ± 0.12 





.77 





,20 


HU 53367 


21. 


04 


± 


0. 


05 


2 


21. 


32 


± 





.30 


Bohlin 


3 


91 «q+0-17 
21.63_ 12 


1 


.78 





.51 


HD 73882 


21. 


.11 


± 


0. 


.08 


1 


21 


.11 


± 





.15 


Ly-a 


4 


91 cn+0.08 

21.59I 007 


1 


.14 





.67 


nn 1 1 n A o 9 

HD 110432 


20 


64 


± 


0. 


04 


1 


20 


85 


± 





15 


Ly-a 


5 


oi on+°- ()8 
21-20_ 007 


1 


.70 





.55 


HD 147888 6 


20 


47 


± 


0. 


05 


2 


21 


.71 


± 





.09 


Ly-a 


2 


21.76 ±0.08 


13 


;.63 





10 


HD 149404 


20. 


79 


± 


0. 


04 


2 


21 


.40 


± 





.14 


Ly-a 


1 


zl '°'-0.10 


1 


.48 





,33 


HD 152236 


20. 


73 


± 


0. 


12 


2 


21 


.77 


± 





.13 


Ly-a 


1 


21.84 ±0.11 


1 


.25 





15 


HD 154368 


21. 


16 


± 


0. 


07 


1 


21. 


.00 


± 





05 


Ly-a 


6 


21.59 ±0.05 


1 


.31 





,74 


HD 167971 


20. 


85 


± 


0. 


12 


1 


21. 


.60 


± 





30 


Ly-a 


7 


21.731°;!* 


2 


.40 





.26 


HD 168076 


20. 


.68 


± 


0. 


08 


1 


21 


.65 


± 





23 


Ly-a 


1 


21.731°;?° 





.84 





.18 


HD 170740 


20. 


86 


± 


0. 


08 


1 


21 


,04 


± 





.15 


Ly-a 


1 


21 41+0-08 


3 


.87 





.57 


HD 179406 


20. 


73 


± 


0. 


07 


2 


21 


.23 


± 





.15 


Shull 


8 


21 44+0.10 


5 


.61 





.39 


HD 185418 


20. 


76 


± 


0. 


05 


1 


21. 


.11 


± 





15 


Ly-a 


4 


91 QQ+0.09 





.83 





.47 


HD 186994 


19. 


59 


± 


0. 


.04 


2 


20 


.90 


± 





.15 


Ly-a 


9 


20.941°;^ 





.11 





.09 


HD 192639 


20. 


69 


± 


0. 


05 


1 


21 


32 


± 





12 


Ly-a 


1 


21.491°;°! 





.90 





.32 


HD 197512 


20. 


.66 


± 


0. 


.05 


1 


21 


.26 


± 





.15 


Ly-a 


4 


21 44+0-11 





.52 





.33 


HD 199579 


20. 


53 


± 


0. 


04 


1 


21 


.04 


± 





.11 


Ly-a 


1 


21.25 ±0.07 





.48 





.38 


HD 203938 


21. 


00 


± 


0. 


06 


1 


21. 


.48 


± 





15 


Ly-a 


4 


21 7n+°' 10 
zi -' u -0.09 


2 


,32 





.40 


HD 206267 


20. 


86 


± 


0. 


.04 


1 


21 


.30 


± 





.15 


Ly-a 


7 


21 54+0-09 


1 


.12 





.42 


HD 207198 


20. 


83 


± 


0. 


04 


1 


21. 


,34 


± 





.17 


Ly-a 


1 


21 59+0- 11 


1 


.15 





.38 


HD 207538 


20. 


.91 


± 


0. 


06 


1 


21 


.34 


± 





.12 


Ly-a 


1 


21 58+0-08 


1 


.40 





.43 


HD 210839 


20. 


84 


± 


0. 


04 


1 


21 


.15 


± 





.10 


Ly-a 


1 


21 45+0-06 


1 


.79 





.49 



'Same as Tabic 2 of Jensen et al. (2005), with inclusion of new lines of sight, minor revisions, and 
correction of typographical errors. 

b Rcfcrences for H 2 : (1) Rachford et al. (2002). (2) B. L. Rachford et al., in preparation. 

c Methods for HI: Ly-a— Profile fitting of the Lyman-a line. Bohlin— N (HI) = 5.8 X 10 21 E B _ V - 2JV(H 2 ) 
from relationship in Bohlin et al. (1978) and further confirmed in Rachford ct al. (2002); errors in 7V(HI) 
assumed to be ±0.30 dex. Shull — ./V(HI) = 5.2 X lO^Bg-v from Shull & Van Stcenburg (1985). 

d References where HI measurements/estimations appear: (1) Diplas & Savage (1994). (2) Cartledge et 
al. (2004). (3) This paper. (4) Fitzpatrick & Massa (1990); errors in 7V(HI) assumed to be ±0.15 dex. (5) 
Rachford ct al. (2001). (6) Snow ct al. (1996). (7) Rachford ct al. (2002). (8) Hanson ct al. (1992). (9) 
Bohlin et al. (1978). 

e See §4.2 for further discussion on these stars. 
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Table 3. Sightlines: Reddening Data 



Sightlinc 


Eb-v (magnitudes) 


Ay (magnitudes) 


7~> 


HD 24534 


0.59 


2.05 


3.47 


HD 27778 


0.37 


1.01 


2.72 


HD 37903 


0.35 


1.28 


3.67 


HD 38087 


0.29 


1.61 


5.57 


HD 40893 


0.46 


1.13 


2.46 


HD 41117 


0.45 


1.23 


2.74 


HD 42087 


0.36 


1.10 


3.06 


HD 46056 


0.50 


1.30 


2.60 


HD 46202 


0.49 


1.39 


2.83 


HD 53367 


0.74 


1.76 


2.38 


HD 73882 


0.70 


2.36 


3.37 


HD 110432 


0.51 


2.02 


3.95 


HD 147888 


0.47 


1.91 


4.06 


HD 149404 


0.68 


2.23 


3.28 


HD 152236 


0.68 


2.24 


3.29 


HD 154368 


0.78 


2.34 


3.00 


HD 167971 


1.08 


3.42 


3.17 


HD 168076 


0.78 


2.77 


3.55 


HD 170740 


0.48 


1.30 


2.71 


HD 179406 


0.33 


0.94 


2.86 


HD 185418 


0.50 


1.16 


2.32 


HD 186994 


0.17 


0.56 


3.29 


HD 192639 


0.66 


1.87 


2.84 


HD 197512 


0.32 


0.75 


2.35 


HD 199579 


0.37 


1.09 


2.95 


HD 203938 


0.74 


2.15 


2.91 


HD 206267 


0.53 


1.41 


2.67 


HD 207198 


0.62 


1.50 


2.42 


HD 207538 


0.64 


1.44 


2.25 


HD 210839 


0.57 


1.58 


2.78 



Table 4. Absorption Lines Observed 



Wavelength (A) Oscillator Strength / Damping Constant 7 Data Used 



951.0792 


1.69 


X 


10- 4 


1.30 


X 


10 8 


FUSE 


951.2948 


2.32 


X 


10- 5 


1.29 


X 


10 8 


FUSE 


953.4152 


1.29 


X 


10- 2 


2.19 


X 


10 8 


FUSE 


953.6549 


2.47 


X 


10- 2 


2.10 


X 


10 8 


FUSE 


959.4937 


5.18 


X 


10- 5 


1.18 


X 


10 8 


FUSE 


1134.1563 


1.46 


X 


10- 2 


1.51 


X 


10 8 


FUSE 


1134.4149 


2.87 


X 


10- 2 


1.49 


X 


10 8 


FUSE 


1134.9803 


4.16 


X 


10~ 2 


1.44 


X 


10 8 


FUSE 


1159.8168 


9.96 


X 


lO" 6 


4.70 


X 


10 8 


FUSE 


1160.9366 


2.75 


X 


10' 6 


4.70 


X 


10 8 


FUSE 


1199.5496 


1.32 


X 


10- 1 


4.07 


X 


10 8 


HST 


1200.2233 


8.69 


X 


ID" 2 


4.02 


X 


10 8 


HST 


1200.7098 


4.32 


X 


10-2 


4.00 


X 


10 8 


HST 
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Table 5. Absorption Lines in the FUSE Wavelength Region Not Used in This Study 



Wavelength (A) Oscillator Strength / Damping Constant 7 Reason Unused a 



952.3034 


2.29 


X 


10- 3 


4.81 


X 


10 7 


H 2 , 01 


952.4148 


1.97 


X 


10- 3 


5.13 


X 


10 7 


H 2 , 01 


952.5227 


5.18 


X 


10- 4 


4.47 


X 


10 7 


H 2 , 01 


953.9699 


3.31 


X 


10- 2 


2.03 


X 


10 8 


H 2 


954.1042 


4.00 


X 


10- 3 


1.41 


X 


10 8 


H 2 


955.2644 


6.92 


X 


10~ 6 


1.13 


X 


10 7 


H 2 


955.8816 


1.91 


X 


10- 5 


1.26 


X 


10 7 


H 2 


960.2014 


1.17 


X 


10- 5 


1.32 


X 


10 8 


Too weak b 


963.9903 


1.24 


X 


io~ 2 


8.55 


X 


10 7 


H 2 


964.6256 


7.90 


X 


10~ 3 


8.31 


X 


10 7 


H 2 


965.0413 


3.86 


X 


10" 3 


8.16 


X 


10 7 


H 2 


1003.3722 


2.81 


X 


10~ 8 


3.99 


X 


10 8 


H 2 , too weak c 


1003.3771 


1.90 


X 


10- 7 


3.99 


X 


10 8 


H 2 , too weak c 



a Describes the reason the line was not observed or conclusively measured, e.g. an 
atomic or H 2 line blends with the NI line or makes the continuum too uncertain to 
measure the NI line, or the NI line is too weak to observe in typical FUSE spectra. 

b The two lines of the 1160 A doublet have weaker oscillator strengths, but the 
S/N is characteristically better redward of 1100 A. The weaker component of a 
doublet at 951 A has an oscillator strength nearly twice that of this line, and is only 
observed in five lines of sight. 

c Observing the 1003 A doublet would require a column density and/or S/N at 
least a combined order of magnitude larger than observed in our best cases. 



Table 6. Measured Equivalent Widths: FUSE Data 



Sightlinc 


W951.1 


W951.3 


W953.4 


W953.7 


W959.5 


^1134.2 


VK1134.4 


W1135.O 


VV1159.S 


160.9 




(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


HD 24534 






68.9 ± 6.5 


83.9 zb 7.0 




99.6 ± 4.4 


121.4 ± 4.4 


136.5 zb 4.4 


11.2 ± 1.4 


5.7 zb 2.1 


HD 27778 












54.3 ± 3.9 


62.5 zb 3.6 


76.2 zb 4.5 


13.2 ± 3.2 


3.3 zb 1.6 


HD 37903 






105.8 ± 11.5 


116.3 zb 11.8 




161.3 zb 2.9 


169.1 ± 3.8 


174.8 zb 3.3 


7.7 ± 1.6 


3.8 zb 3.2 


HD 38087 












173. 1 zb 6.9 


183.7 zb 6.6 


189.6 zb 7.1 


7.8 ± 4.0 




HD 40893 






260.8 ± 88.9 


294.1 zb 126.1 


43.0 zb 12.7 


293.1 zb 13.6 


316.3 ± 24.1 


332.9 zb 31.1 


28.8 =t 4.3 


4.4 zb 3.1 


HD 41117 












146.7 ± 11.5 


170.3 zb 5. 1 


173.4 zb 17.6 


32.8 ± 19.5 




HD 42087 






130.1 ± 33.7 


100.7 ± 28.8 


54.1 zb 16.4 


148.7 ± 3.5 


180.2 zb 5.5 


250.8 ± 9.2 


17.2 zb 3.6 




HD 46056 












152.9 zb 11.2 


165.9 ± 11.3 


172.8 ± 12.0 


25.3 zb 9.9 




HD 46202 


57.8 ± 18.8 




106.7 zb 12.0 


118.9 ± 13.8 


65.7 ± 22.3 


138.1 zb 14.4 


157.2 ± 13.7 


166.8 ± 14.2 




7.8 zb 3.5 


HD 53367 












133.7 ± 19.2 


159.0 ± 17.3 


178.9 ± 15.0 


24.4 zb 4.3 




HD 73882 












191.7 zb 6.7 


210.0 ± 6.8 


219.1 ± 6.3 


22.6 zb 1.8 


8.2 zb 2.1 


HD 110432 


30.0 ± 3.8 


7.2 ± 2.4 


80.5 zb 3.9 


92.3 ± 5.7 


23.2 zb 4.0 


105.5 zb 1.0 


118.1 ±1.1 


123.9 ± 1.1 


14.3 zb 1.5 


4.3 zb 1.1 


HD 147888 






84.0 ± 22.9 


78.7 ± 26.5 


49.6 ± 20.4 


116.5 zb 5.0 


132.9 zb 5.4 


124.8 ± 13.2 


16.5 zb 3.7 


5.8 ± 3.5 


HD 149404 


62.8 ± 5.0 


27.7 dz 4.3 


131.7 zb 8.7 


140.6 ± 9.2 


42.9 ± 4.2 


171.1 ± 6.4 


179.7 zb 6.5 


184.5 ± 6.5 


34.1 zb 1.0 


11.2 zb 1.7 


HD 152236 












296.0 ± 27.3 


321.9 ± 43.1 


340.6 ± 49.5 


20.4 zb 3.5 




HD 154368 












195.8 ± 14.2 


206.3 ± 16.8 


213.2 ± 19.8 


41.7 zb 7.8 




HD 167971 












223.8 ± 23.6 


236.1 ± 23.8 


242.7 ± 23.9 


55.9 zb 13.0 




HD 168076 












250.3 ± 19.0 


277.7 ± 17.2 


291.8 ± 17.1 


38.7 zb 8.9 




HD 170740 






58.5 ± 34.9 


85.2 ± 24.4 


66.1 zb 21.3 


97.8 ± 3.7 


126.5 ± 4.4 


123.0 ± 4.5 




10.0 =b 4.0 


HD 179406 












114.1 ± 5.6 


121.6 ± 5.4 


132.1 ± 5.7 


16.6 zb 10.6 




HD 185418 






79.3 ± 12.7 


95.0 ± 14.7 




107.6 ± 1.1 


122.6 ± 1.3 


133.4 ± 1.4 


15.8 zb 2.2 


2.6 zb 2.1 


HD 186994 


22.1 ± 11.5 




224.8 ± 12.2 


243.2 ± 7.3 




277.0 ± 9.5 


289.8 ± 10.2 


297.0 ± 10.7 






HD 192639 


48.7 ± 14.1 


38.0 ± 15.0 


159.3 ± 13.3 


178.4 ± 14.6 


45.1 zb 14.7 


246.4 ± 7.0 


230.7 ± 5.0 


251.0 ± 6.8 


16.6 zb 3.6 


4.8 zb 3.4 


HD 197512 






142.3 ± 16.8 


166.2 ± 10.4 


76.5 zb 17.8 


124.1 ± 5.6 


139.7 ± 6.2 


137.0 ± 9.3 


6.4 =b 1.9 




HD 199579 


53.3 ± 3.8 


15.6 ± 3.5 


159.3 ± 13.3 


178.4 ± 14.6 


25.0 ± 3.3 


179.6 ± 2.0 


203.8 ± 2.6 


210.2 ± 3.1 


19.7 zb 1.4 


6.6 zb 2.2 


HD 203938 












131.4 ± 12.3 


139.9 ± 12.1 


144.6 ± 12.4 


29.1 zb 7.8 




HD 206267 


80.0 ± 42.4 




100.9 ± 18.4 


129.1 ± 10.7 


26.0 zb 14.4 


156.3 ± 2.4 


179.9 ± 1.6 


159.0 ± 6.4 


15.3 zb 2.6 


6.3 zb 2.4 


HD 207198 


62.7 ± 24.2 




126.6 ± 33.5 


108.1 ± 17.1 


63.4 =b 21.0 


135.1 ± 5.9 


142.6 ± 5.8 


146.8 ± 6.0 


22.7 zb 2.6 


5.4 =b 2.0 


HD 207538 






110.9 ± 44.4 


125.6 ± 30.4 




125.3 ± 6.4 


136.3 ± 6.5 


141.9 ± 6.2 


27.7 zb 3.6 




HD 210839 


55.3 ± 4.8 


20.7 ± 4.6 


133.2 ± 8.6 


162.1 ± 6.5 


22.4 ± 3.8 


186.1 ± 2.5 


193.2 ± 3.0 


200.0 ± 3.8 


18.6 zb 1.1 


4.6 zb 0.9 



Table 7. Measured Equivalent Widths: HST Data 



Sightline W1199.5 ^1200.2 W1200.7 

(mA) (mA) (mA) 



HD 27778 
HD 37903 
HD 147888 
HD 185418 
HD 192639 
HD 206267 
HD 207198 
HD 207538 
HD 210839 



110.7 ±6.4 

224.1 ± 7.3 
195.6 ± 11.3 

180.2 ±5.5 
323.5 ±68.2 

216.8 ±28.2 
207.1 ±6.7 

182.9 ±8.6 
251.9 ±20.9 



103.5 ±7.1 

215.6 ±8.0 
182.5 ±11.4 
170.4 ±6.1 

303.3 ±49.2 

201.4 ±20.6 
198.9 ±6.9 

172.7 ±8.4 
235.0 ± 15.4 



89.8 ±7.5 

200.8 ± 8.4 
157.3 ± 12.5 

152.3 ±6.3 
277.7 ±25.4 

182.9 ± 11.0 
184.6 ±6.6 
153.0 ±8.9 
214.0 ±8.1 
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Table 8. Column Density Results 



Sightrlinc 


log NI 


6-V3,lllC 


log NH 


N/H (nnml 


Method 3, 


HD 24534 


17 99 +0.09 
1 1 - zz _o.08 


4 9+0-9 


-4.12 + 0.09 




Uniform Error COG b 


HD 27778 


17.11 ± 0.16 


2.7 + 0.4 


-4 29+ 16 
y — 18 


5llf 7 


COG 


HD 37903 


16.89±g-J| 


8.8 + 0.8 


__ 4 fil +0.38 


24lit 


COG 


HD 38087 


16.8312$ 


9.6 + 0.7 


^•^ u -0.49 


40lf 6 


COG 


HD 40893 


17.38l°;?i 




-4 91+0- 22 
4 - zl -0.23 


6111 


COG 


HD 41117 


18.53l?-y 


4.6l 3 2 2 


-s m+o 18 

o.Ul_ 2 99 


976lf 7 4 


COG 


HD 42087 


17.23±g;i* 


8-3l?i 


-4 26+ ' 15 
^• zu -0.18 


54+23 


Uniform Error COG 


HD 46056 


17 40+0.27 


7.5 + 0.6 


_ 4 ,,+0.28 
4 - ii -0.30 


78+69 
'8-38 


COG 


HD 46202 


17 4 2+0.46 


6.9 ± 1.2 


-4 26+ ' 46 
4 - ZD -0.47 


55+104 
00-36 


COG 


HD 53367 


17 4Q+0.07 


7.2 + 0.5 


-4 93+0-12 
4 - ZJ -0.31 


581^ 


COG b 


HD 73882 


17.36 ± 0.11 


Q Q+0.6 


-4 23+ ' 12 


591J? 


COG 


HD 110432 


17.09l°;ii 


5.5 + 0.3 


— 4 1 1 + - 19 


77 +44 
"-30 


COG 


HD 147888 


17.38l°;f7 




-4.38l2-jt 


41+50 


COG 


HD 149404 


17.561^1 


7.8 + 0.6 


-4.oii2;i° 


9611 


COG 


HD 152236 


17.27 + 0.08 


15.9 ± 1.4 


-4.57±8;i| 


26+ 9 


COG 


HD 154368 


17.67l° ;i? 


9 2+°' 7 


-3 O9+0- 14 
• 3 - yz -0.13 




COG 


HD 167971 


ir.83±g;ig 


10 3+ 11 


-3 Q0+0- 21 
• 3 - yu -0.69 


125199 


COG 


HD 168076 




12.7 + 0.8 


_4 u+0- 17 


71 + 34 
'1_44 


COG 


HD 170740 


17.65l°;l 6 7 


4.8l°J 


o 7f?+0.56 
— a -' D -0.41 


17 5+462 
1,o -107 


COG 


HD 179406 


17 97 +i.oo 


5.6+ij 


_ 4 17 +1.00 
-0.89 


67l 6 5 ° 4 


COG 


HD 185418 


17.04l° ;g 


5-91?;° 


_ 4 Qr+0.25 
1.OO_ 27 


44+H 


Uniform Error COG 


HD 186994 


16.59lJ;° 3 2 


i6.8l 2 3 ;! 


_ 4 35+103 
^•" ; ' o -0.68 


44+ 4 f 


COG 


HD 192639 




12.0 ± 1.4 


_4 31+0.48 
4 - Ji -0.44 


49+2? 


COG 


HD 197512 


16.9612:19 


7.3 + 0.6 


4 4S+0.22 
-4.48_ 26 


33±?1 


COG 


HD 199579 


17.16 ± 0.10 


o 7+0-4 
a -'-0.3 


-4 nq+o- 11 

4 - ua -0.14 


811^ 


COG 


HD 203938 


17.5312;1 9 7 


6.1 + 0.6 


_ 4 17+O.20 
4 - 1( -0.23 


67lf 7 


COG 


HD 206267 


17.161^ 


8.5 + 0.6 


4 OO+0.27 
-4.38_ 34 


4111 


COG 


HD 207198 


1 ' -0.19 


6.7 + 0.7 


_ 4 14+0.19 
4 - 14 -0.27 


72l|2 


COG 


HD 207538 


17.39 ± 0.17 


6.0 + 0.6 


_4 19+O.I8 




COG 


HD 210839 


I7.i9l° ;« 


9.5 + 0.5 


-4.26 + 0.14 


55+20 
00 -15 


COG 



a Methods: COG — Standard curve-of-growth method. Uniform Error COG — Curve-of-growth 
method with fractional errors averaged. 

b This solution is taken from a local minimum in the \ 2 array (from the stated curve-of-growth 
method) that is not the global minimum; the global minimum is ruled out based on inconsistency 
with the line profiles. See §4.3 for further discussion. 



